InĀ [1]:
# Import libraries
import os
import pathlib

import hvplot.pandas
import hvplot.xarray
import holoviews as hv

import geopandas as gpd
import rioxarray as rxr
import xarray as xr
import pandas as pd

from pathlib import Path
from glob import glob
import earthpy
import earthpy.api.appeears as eaapp
InĀ [2]:
# Make project folder
project = earthpy.Project(
    'Gila River Vegetation', dirname='Gila_Vegetation_Data')
project.get_data()

# Load in the boundary data
aitsn_gdf = gpd.read_file(project.project_dir / 'tl_2020_us_aitsn')
InĀ [3]:
# Select and merge the subdivisions you want
prir_gdf = aitsn_gdf[aitsn_gdf["AIANNHCE"]== '2810']
prir_gdf #for Ogala Lakota Nation

# Plot the results with web tile images
prir_gdf.plot()
Out[3]:
<Axes: >
No description has been provided for this image
InĀ [4]:
# Select and merge the subdivisions you want
prir_gdf = prir_gdf.loc[prir_gdf.AIANNHCE=='2810'].dissolve()
# Plot the results with web tile images
prir_gdf.hvplot(
    geo=True, tiles='EsriImagery',
    fill_color=None, line_color='black',
    title="Pine Ridge Indian Community",
    frame_width=500)
Out[4]:
InĀ [5]:
# Initialize AppeearsDownloader for MODIS NDVI data
ndvi_downloader = eaapp.AppeearsDownloader(
    download_key='pric-ndvi',
    project=project,
    product='MOD13Q1.061',
    layer='_250m_16_days_NDVI',
    start_date="06-01",
    end_date="09-01",
    recurring=True,
    year_range=[2016, 2022],
    polygon=prir_gdf
)
InĀ [6]:
# Download the prepared download -- this can take some time!
ndvi_downloader.download_files(cache=True)
Out[6]:
<generator object Path.rglob at 0x74591c2b3670>
InĀ [7]:
# Get a list of NDVI tif file paths
project.project_dir = '/workspaces/data/gila_vegetation_data/pric-ndvi'
project.project_dir

ndvi_paths = sorted(glob(os.path.join(project.project_dir, '**', '*NDVI*.tif')))
ndvi_paths
Out[7]:
['/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2016241000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2017241000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2018241000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2019241000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2020241000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2021241000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022145000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022161000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022177000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022193000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022209000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022225000000_aid0001.tif',
 '/workspaces/data/gila_vegetation_data/pric-ndvi/MOD13Q1.061_2016138_to_2022244/MOD13Q1.061__250m_16_days_NDVI_doy2022241000000_aid0001.tif']
InĀ [8]:
doy_start = -25
doy_end = -18

# Loop through each NDVI image
ndvi_das = []
for ndvi_path in ndvi_paths:
    ndvi_path = Path(ndvi_path)

    # Get date from file name
    doy = ndvi_path.name[doy_start:doy_end]
    date = pd.to_datetime(doy, format='%Y%j')

    # Open dataset
    da = rxr.open_rasterio(ndvi_path, mask_and_scale=True).squeeze()

    # Add date dimension and clean up metadata
    da = da.assign_coords({'date': date})
    da = da.expand_dims({'date': 1})
    da.name = 'NDVI'

    # Prepare for concatenation
    ndvi_das.append(da)

ndvi_das
Out[8]:
[<xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3466    , 0.34719998, 0.34329998, ..., 0.7053    ,
          0.4621    , 0.5223    ],
         [0.3511    , 0.3241    , 0.3265    , ..., 0.546     ,
          0.4648    , 0.4648    ],
         [0.3463    , 0.3637    , 0.3637    , ..., 0.5022    ,
          0.4648    , 0.4648    ],
         ...,
         [0.498     , 0.54609996, 0.55869997, ..., 0.6034    ,
          0.6034    , 0.4711    ],
         [0.5054    , 0.4606    , 0.5184    , ..., 0.68509996,
          0.61869997, 0.5227    ],
         [0.46789998, 0.5051    , 0.5342    , ..., 0.67039996,
          0.5625    , 0.5503    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-05-24
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.36389998, 0.3691    , 0.3685    , ..., 0.5238    ,
          0.4243    , 0.40219998],
         [0.36589998, 0.36589998, 0.3522    , ..., 0.4353    ,
          0.4243    , 0.4243    ],
         [0.361     , 0.3703    , 0.3703    , ..., 0.4353    ,
          0.4622    , 0.4622    ],
         ...,
         [0.3988    , 0.44079998, 0.4702    , ..., 0.7672    ,
          0.7672    , 0.5457    ],
         [0.41149998, 0.40719998, 0.504     , ..., 0.6691    ,
          0.6019    , 0.57879996],
         [0.4252    , 0.3804    , 0.43899998, ..., 0.6691    ,
          0.55469996, 0.5628    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-06-09
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.31039998, 0.29729998, 0.3117    , ..., 0.38529998,
          0.384     , 0.4831    ],
         [0.2959    , 0.29099998, 0.29099998, ..., 0.3557    ,
          0.4274    , 0.4274    ],
         [0.2895    , 0.2933    , 0.2933    , ..., 0.4597    ,
          0.4143    , 0.4143    ],
         ...,
         [0.4112    , 0.42279997, 0.415     , ..., 0.68      ,
          0.68      , 0.64419997],
         [0.4141    , 0.4141    , 0.4244    , ..., 0.701     ,
          0.6992    , 0.6871    ],
         [0.43789998, 0.38529998, 0.4244    , ..., 0.7391    ,
          0.6992    , 0.6871    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-06-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3093    , 0.3217    , 0.3324    , ..., 0.49539998,
          0.3825    , 0.4233    ],
         [0.3093    , 0.3024    , 0.30609998, ..., 0.3964    ,
          0.4278    , 0.4278    ],
         [0.29909998, 0.3024    , 0.3024    , ..., 0.40899998,
          0.4276    , 0.4276    ],
         ...,
         [0.3622    , 0.36879998, 0.36949998, ..., 0.6454    ,
          0.6454    , 0.6723    ],
         [0.3583    , 0.3682    , 0.3958    , ..., 0.62619996,
          0.658     , 0.65309995],
         [0.37669998, 0.3487    , 0.391     , ..., 0.70199996,
          0.6875    , 0.6454    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-07-11
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3283    , 0.33269998, 0.3604    , ..., 0.5243    ,
          0.35439998, 0.44309998],
         [0.3285    , 0.33269998, 0.332     , ..., 0.392     ,
          0.4069    , 0.4069    ],
         [0.3285    , 0.33429998, 0.33429998, ..., 0.4109    ,
          0.40379998, 0.40379998],
         ...,
         [0.38509998, 0.40289998, 0.39909998, ..., 0.7237    ,
          0.7237    , 0.7237    ],
         [0.38509998, 0.40289998, 0.406     , ..., 0.7371    ,
          0.7435    , 0.6936    ],
         [0.38819999, 0.38819999, 0.4192    , ..., 0.7215    ,
          0.6932    , 0.6971    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-07-27
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.4312    , 0.4454    , 0.46269998, ..., 0.51629996,
          0.40039998, 0.46069998],
         [0.4497    , 0.46429998, 0.48499998, ..., 0.4511    ,
          0.48499998, 0.48499998],
         [0.4782    , 0.46429998, 0.46429998, ..., 0.4511    ,
          0.4421    , 0.4421    ],
         ...,
         [0.3585    , 0.3881    , 0.40039998, ..., 0.75699997,
          0.75699997, 0.6455    ],
         [0.3874    , 0.3768    , 0.40039998, ..., 0.7445    ,
          0.7624    , 0.6809    ],
         [0.39389998, 0.3768    , 0.3906    , ..., 0.6595    ,
          0.66069996, 0.60969996]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-08-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.4518    , 0.46109998, 0.4966    , ..., 0.6502    ,
          0.70089996, 0.4815    ],
         [0.49609998, 0.5175    , 0.5164    , ..., 0.6502    ,
          0.48049998, 0.48049998],
         [0.5043    , 0.5773    , 0.5773    , ..., 0.5153    ,
          0.5139    , 0.5139    ],
         ...,
         [0.3437    , 0.3699    , 0.3838    , ..., 0.7194    ,
          0.7194    , 0.6369    ],
         [0.34759998, 0.3489    , 0.37399998, ..., 0.6602    ,
          0.5949    , 0.6369    ],
         [0.3533    , 0.3613    , 0.3868    , ..., 0.6469    ,
          0.5949    , 0.5699    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2016-08-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.5789    , 0.5726    , 0.5541    , ..., 0.6509    ,
          0.6448    , 0.42049998],
         [0.5986    , 0.5726    , 0.56799996, ..., 0.5211    ,
          0.41869998, 0.41869998],
         [0.5986    , 0.5986    , 0.5986    , ..., 0.4094    ,
          0.4671    , 0.4671    ],
         ...,
         [0.5906    , 0.589     , 0.69159997, ..., 0.5413    ,
          0.5413    , 0.4314    ],
         [0.59999996, 0.589     , 0.64559996, ..., 0.5855    ,
          0.5254    , 0.4942    ],
         [0.529     , 0.4955    , 0.57049996, ..., 0.6074    ,
          0.5855    , 0.5513    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-05-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.5826    , 0.58129996, 0.56299996, ..., 0.5846    ,
          0.5056    , 0.5056    ],
         [0.5894    , 0.5755    , 0.5755    , ..., 0.5222    ,
          0.4411    , 0.4411    ],
         [0.6045    , 0.6209    , 0.6209    , ..., 0.45749998,
          0.44129997, 0.44129997],
         ...,
         [0.5087    , 0.5392    , 0.5909    , ..., 0.6357    ,
          0.6357    , 0.5316    ],
         [0.5837    , 0.5281    , 0.5909    , ..., 0.6655    ,
          0.6357    , 0.539     ],
         [0.5559    , 0.5281    , 0.5599    , ..., 0.685     ,
          0.6096    , 0.6096    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-06-10
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3909    , 0.38549998, 0.38329998, ..., 0.4516    ,
          0.3581    , 0.4226    ],
         [0.4014    , 0.41459998, 0.41009998, ..., 0.33359998,
          0.3839    , 0.3839    ],
         [0.4382    , 0.43809998, 0.43809998, ..., 0.3931    ,
          0.39569998, 0.39569998],
         ...,
         [0.433     , 0.4626    , 0.4626    , ..., 0.5652    ,
          0.5652    , 0.4622    ],
         [0.433     , 0.4278    , 0.45189998, ..., 0.57229996,
          0.5378    , 0.4964    ],
         [0.4219    , 0.4219    , 0.4227    , ..., 0.57049996,
          0.5459    , 0.5378    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-06-26
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.305     , 0.3032    , 0.3141    , ..., 0.31149998,
          0.3116    , 0.3191    ],
         [0.31579998, 0.31579998, 0.3141    , ..., 0.3105    ,
          0.3249    , 0.3249    ],
         [0.3215    , 0.327     , 0.327     , ..., 0.32889998,
          0.343     , 0.343     ],
         ...,
         [0.34039998, 0.3338    , 0.3552    , ..., 0.47739998,
          0.47739998, 0.3811    ],
         [0.3322    , 0.3338    , 0.3552    , ..., 0.5191    ,
          0.4405    , 0.42479998],
         [0.3463    , 0.3322    , 0.3353    , ..., 0.4861    ,
          0.44399998, 0.4554    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-07-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3005    , 0.30249998, 0.30249998, ..., 0.3351    ,
          0.325     , 0.3716    ],
         [0.2855    , 0.2947    , 0.30249998, ..., 0.31329998,
          0.31329998, 0.31329998],
         [0.3129    , 0.3058    , 0.3058    , ..., 0.3006    ,
          0.34309998, 0.34309998],
         ...,
         [0.3198    , 0.34309998, 0.3193    , ..., 0.5688    ,
          0.5688    , 0.4678    ],
         [0.3198    , 0.31      , 0.3256    , ..., 0.52      ,
          0.47309998, 0.47309998],
         [0.2891    , 0.3071    , 0.3213    , ..., 0.4605    ,
          0.4623    , 0.47309998]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-07-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.2905    , 0.2924    , 0.29099998, ..., 0.3283    ,
          0.327     , 0.2999    ],
         [0.3012    , 0.30719998, 0.305     , ..., 0.35819998,
          0.331     , 0.331     ],
         [0.2996    , 0.3145    , 0.3145    , ..., 0.35819998,
          0.3464    , 0.3464    ],
         ...,
         [0.2881    , 0.27539998, 0.2733    , ..., 0.7382    ,
          0.7382    , 0.6105    ],
         [0.27809998, 0.2571    , 0.2733    , ..., 0.6545    ,
          0.589     , 0.5243    ],
         [0.2963    , 0.2543    , 0.2733    , ..., 0.6261    ,
          0.58599997, 0.5243    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-08-13
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.2607    , 0.2697    , 0.2814    , ..., 0.39139998,
          0.2777    , 0.2943    ],
         [0.2645    , 0.2697    , 0.27699998, ..., 0.2951    ,
          0.326     , 0.326     ],
         [0.2837    , 0.2848    , 0.2848    , ..., 0.3057    ,
          0.3293    , 0.3293    ],
         ...,
         [0.28939998, 0.28939998, 0.28939998, ..., 0.7254    ,
          0.7254    , 0.58239996],
         [0.2825    , 0.2825    , 0.2825    , ..., 0.6572    ,
          0.6685    , 0.6325    ],
         [0.2853    , 0.2853    , 0.2853    , ..., 0.6511    ,
          0.6181    , 0.5821    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2017-08-29
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.6132    , 0.6095    , 0.6142    , ..., 0.7292    ,
          0.6143    , 0.44979998],
         [0.61759996, 0.6277    , 0.62689996, ..., 0.59529996,
          0.59529996, 0.59529996],
         [0.6306    , 0.67689997, 0.67689997, ..., 0.5908    ,
          0.4622    , 0.4622    ],
         ...,
         [0.4337    , 0.44919997, 0.54609996, ..., 0.72639996,
          0.72639996, 0.5654    ],
         [0.4423    , 0.44919997, 0.5106    , ..., 0.655     ,
          0.655     , 0.5855    ],
         [0.44579998, 0.4427    , 0.4583    , ..., 0.6666    ,
          0.6497    , 0.5669    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-05-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.6468    , 0.6202    , 0.6062    , ..., 0.5879    ,
          0.5044    , 0.39909998],
         [0.6472    , 0.6472    , 0.6486    , ..., 0.598     ,
          0.5044    , 0.5044    ],
         [0.64059997, 0.6472    , 0.6472    , ..., 0.5793    ,
          0.47      , 0.47      ],
         ...,
         [0.511     , 0.5093    , 0.58239996, ..., 0.7       ,
          0.7       , 0.6074    ],
         [0.511     , 0.5093    , 0.5744    , ..., 0.7119    ,
          0.69549996, 0.6753    ],
         [0.5033    , 0.5033    , 0.55439997, ..., 0.7342    ,
          0.7342    , 0.6206    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-06-10
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.7111    , 0.68979996, 0.69299996, ..., 0.6189    ,
          0.47309998, 0.5664    ],
         [0.68479997, 0.704     , 0.69299996, ..., 0.58169997,
          0.5606    , 0.5606    ],
         [0.6835    , 0.704     , 0.704     , ..., 0.5893    ,
          0.5676    , 0.5676    ],
         ...,
         [0.5001    , 0.5451    , 0.5751    , ..., 0.7869    ,
          0.7869    , 0.6936    ],
         [0.5001    , 0.4989    , 0.5502    , ..., 0.73649997,
          0.7468    , 0.74189997],
         [0.5115    , 0.5169    , 0.5974    , ..., 0.7795    ,
          0.8081    , 0.76129997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-06-26
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.7389    , 0.70339996, 0.67649996, ..., 0.3585    ,
          0.3847    , 0.457     ],
         [0.7389    , 0.7242    , 0.7213    , ..., 0.4709    ,
          0.42589998, 0.42589998],
         [0.7312    , 0.7076    , 0.7076    , ..., 0.4766    ,
          0.42589998, 0.42589998],
         ...,
         [0.4865    , 0.4883    , 0.4956    , ..., 0.7786    ,
          0.7786    , 0.7309    ],
         [0.4574    , 0.4696    , 0.511     , ..., 0.7382    ,
          0.7481    , 0.7206    ],
         [0.47149998, 0.4513    , 0.5002    , ..., 0.766     ,
          0.78849995, 0.76489997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-07-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.72569996, 0.719     , 0.7064    , ..., 0.3296    ,
          0.28309998, 0.343     ],
         [0.7283    , 0.7295    , 0.71559995, ..., 0.3534    ,
          0.3618    , 0.3618    ],
         [0.723     , 0.7247    , 0.7247    , ..., 0.37359998,
          0.371     , 0.371     ],
         ...,
         [0.5       , 0.52349997, 0.5018    , ..., 0.70409995,
          0.70409995, 0.70849997],
         [0.5       , 0.49229997, 0.49139997, ..., 0.7436    ,
          0.7436    , 0.73859996],
         [0.5417    , 0.5053    , 0.5178    , ..., 0.80359995,
          0.7733    , 0.73859996]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-07-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.509     , 0.513     , 0.49929997, ..., 0.3931    ,
          0.3401    , 0.3621    ],
         [0.515     , 0.4919    , 0.482     , ..., 0.40609998,
          0.3845    , 0.3845    ],
         [0.52599996, 0.5155    , 0.5155    , ..., 0.41259998,
          0.3937    , 0.3937    ],
         ...,
         [0.4631    , 0.4752    , 0.4754    , ..., 0.70239997,
          0.70239997, 0.6732    ],
         [0.4631    , 0.46679997, 0.48209998, ..., 0.7499    ,
          0.69659996, 0.6732    ],
         [0.41529998, 0.4748    , 0.48479998, ..., 0.7098    ,
          0.6987    , 0.6732    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-08-13
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.4355    , 0.44779998, 0.43039998, ..., 0.59779996,
          0.4434    , 0.48      ],
         [0.4549    , 0.4369    , 0.36699998, ..., 0.5229    ,
          0.48299998, 0.48299998],
         [0.4549    , 0.4253    , 0.4253    , ..., 0.42459998,
          0.48299998, 0.48299998],
         ...,
         [0.31129998, 0.3681    , 0.37579998, ..., 0.671     ,
          0.671     , 0.4966    ],
         [0.3254    , 0.31739998, 0.37579998, ..., 0.6243    ,
          0.671     , 0.5768    ],
         [0.36139998, 0.35529998, 0.3737    , ..., 0.6243    ,
          0.6026    , 0.5692    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2018-08-29
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.65819997, 0.6678    , 0.6892    , ..., 0.57879996,
          0.75      , 0.33589998],
         [0.7073    , 0.728     , 0.7294    , ..., 0.57879996,
          0.3129    , 0.3129    ],
         [0.73179996, 0.7689    , 0.7689    , ..., 0.2753    ,
          0.26299998, 0.26299998],
         ...,
         [0.608     , 0.6277    , 0.6519    , ..., 0.65239996,
          0.65239996, 0.5064    ],
         [0.587     , 0.57519996, 0.62939996, ..., 0.6583    ,
          0.559     , 0.47059998],
         [0.60499996, 0.6742    , 0.7101    , ..., 0.5299    ,
          0.43789998, 0.4781    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-05-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.7877    , 0.7877    , 0.796     , ..., 0.6807    ,
          0.7999    , 0.42229998],
         [0.7877    , 0.7877    , 0.796     , ..., 0.6807    ,
          0.4692    , 0.4692    ],
         [0.76629996, 0.81009996, 0.81009996, ..., 0.40989998,
          0.3741    , 0.3741    ],
         ...,
         [0.648     , 0.709     , 0.69729996, ..., 0.7769    ,
          0.7769    , 0.6732    ],
         [0.6757    , 0.63659996, 0.70769995, ..., 0.70769995,
          0.6724    , 0.6274    ],
         [0.6674    , 0.68299997, 0.74909997, ..., 0.7446    ,
          0.75      , 0.6752    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-06-10
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.7758    , 0.7407    , 0.7525    , ..., 0.77629995,
          0.77629995, 0.4684    ],
         [0.7758    , 0.77809995, 0.77809995, ..., 0.6437    ,
          0.5998    , 0.5998    ],
         [0.75869995, 0.7974    , 0.7974    , ..., 0.6437    ,
          0.5998    , 0.5998    ],
         ...,
         [0.63339996, 0.71919996, 0.7129    , ..., 0.6836    ,
          0.6836    , 0.6597    ],
         [0.66569996, 0.6307    , 0.69409996, ..., 0.7341    ,
          0.79459995, 0.7116    ],
         [0.66029996, 0.6641    , 0.72099996, ..., 0.80249995,
          0.79459995, 0.7227    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-06-26
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.72789997, 0.719     , 0.7346    , ..., 0.7734    ,
          0.7298    , 0.611     ],
         [0.73719996, 0.73249996, 0.7301    , ..., 0.6685    ,
          0.5451    , 0.5451    ],
         [0.7044    , 0.73249996, 0.73249996, ..., 0.6116    ,
          0.53459996, 0.53459996],
         ...,
         [0.58599997, 0.5979    , 0.6464    , ..., 0.7653    ,
          0.7653    , 0.68299997],
         [0.62369996, 0.5979    , 0.5979    , ..., 0.6641    ,
          0.7132    , 0.6724    ],
         [0.62369996, 0.5826    , 0.6562    , ..., 0.7183    ,
          0.766     , 0.7241    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-07-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.6838    , 0.6838    , 0.67829996, ..., 0.5819    ,
          0.5829    , 0.4737    ],
         [0.6796    , 0.6691    , 0.67829996, ..., 0.7187    ,
          0.5829    , 0.5829    ],
         [0.6796    , 0.6691    , 0.6691    , ..., 0.7187    ,
          0.6087    , 0.6087    ],
         ...,
         [0.5869    , 0.5869    , 0.6155    , ..., 0.70739996,
          0.70739996, 0.6824    ],
         [0.615     , 0.5775    , 0.6155    , ..., 0.709     ,
          0.7181    , 0.7181    ],
         [0.63629997, 0.5775    , 0.65349996, ..., 0.7604    ,
          0.7574    , 0.67069995]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-07-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.62259996, 0.62189996, 0.6195    , ..., 0.7321    ,
          0.5875    , 0.57629997],
         [0.61869997, 0.61869997, 0.5937    , ..., 0.73139995,
          0.6785    , 0.6785    ],
         [0.6246    , 0.61939996, 0.61939996, ..., 0.73139995,
          0.66819996, 0.66819996],
         ...,
         [0.5445    , 0.5433    , 0.57269996, ..., 0.76269996,
          0.76269996, 0.7407    ],
         [0.5186    , 0.5186    , 0.5639    , ..., 0.76089996,
          0.756     , 0.7407    ],
         [0.52309996, 0.52309996, 0.54109997, ..., 0.7804    ,
          0.756     , 0.6846    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-08-13
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.57049996, 0.5568    , 0.5457    , ..., 0.6229    ,
          0.5591    , 0.5591    ],
         [0.5546    , 0.5148    , 0.5457    , ..., 0.6552    ,
          0.5667    , 0.5667    ],
         [0.5609    , 0.5527    , 0.5527    , ..., 0.6321    ,
          0.6005    , 0.6005    ],
         ...,
         [0.36089998, 0.3884    , 0.40269998, ..., 0.6951    ,
          0.6951    , 0.5972    ],
         [0.3798    , 0.3724    , 0.38279998, ..., 0.643     ,
          0.643     , 0.6217    ],
         [0.4059    , 0.3703    , 0.41329998, ..., 0.70419997,
          0.652     , 0.5803    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2019-08-29
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3913    , 0.4057    , 0.4459    , ..., 0.807     ,
          0.5846    , 0.5846    ],
         [0.4148    , 0.4057    , 0.4263    , ..., 0.69769996,
          0.5744    , 0.5744    ],
         [0.4148    , 0.4154    , 0.4154    , ..., 0.5085    ,
          0.4807    , 0.4807    ],
         ...,
         [0.4254    , 0.4364    , 0.45139998, ..., 0.5872    ,
          0.5872    , 0.522     ],
         [0.4247    , 0.4247    , 0.44759998, ..., 0.65169996,
          0.60069996, 0.522     ],
         [0.4385    , 0.4247    , 0.44759998, ..., 0.5708    ,
          0.5304    , 0.523     ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-05-24
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.406     , 0.4026    , 0.4113    , ..., 0.7705    ,
          0.6789    , 0.53639996],
         [0.38869998, 0.4008    , 0.40449998, ..., 0.6033    ,
          0.5376    , 0.5376    ],
         [0.41529998, 0.4095    , 0.4095    , ..., 0.5599    ,
          0.5376    , 0.5376    ],
         ...,
         [0.46449998, 0.4664    , 0.4684    , ..., 0.704     ,
          0.704     , 0.539     ],
         [0.43289998, 0.43289998, 0.4491    , ..., 0.722     ,
          0.6684    , 0.6313    ],
         [0.42929998, 0.4428    , 0.4428    , ..., 0.69589996,
          0.6738    , 0.61399996]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-06-09
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.38259998, 0.3707    , 0.3775    , ..., 0.6134    ,
          0.4763    , 0.4763    ],
         [0.3672    , 0.3682    , 0.3642    , ..., 0.60389996,
          0.5252    , 0.5252    ],
         [0.3867    , 0.3915    , 0.3915    , ..., 0.6116    ,
          0.5844    , 0.5844    ],
         ...,
         [0.3617    , 0.39299998, 0.4111    , ..., 0.66569996,
          0.66569996, 0.5837    ],
         [0.39499998, 0.3753    , 0.39389998, ..., 0.6979    ,
          0.68009996, 0.64599997],
         [0.377     , 0.3753    , 0.39389998, ..., 0.6817    ,
          0.6789    , 0.64699996]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-06-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3636    , 0.3665    , 0.371     , ..., 0.61649996,
          0.61649996, 0.4881    ],
         [0.3545    , 0.3545    , 0.3464    , ..., 0.687     ,
          0.5905    , 0.5905    ],
         [0.3656    , 0.3656    , 0.3656    , ..., 0.687     ,
          0.5905    , 0.5905    ],
         ...,
         [0.3112    , 0.3505    , 0.3505    , ..., 0.7542    ,
          0.7542    , 0.59419996],
         [0.3112    , 0.332     , 0.3335    , ..., 0.6685    ,
          0.71459997, 0.6817    ],
         [0.3278    , 0.3278    , 0.3193    , ..., 0.66249996,
          0.69      , 0.6885    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-07-11
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3491    , 0.3355    , 0.32889998, ..., 0.6778    ,
          0.67359996, 0.49089998],
         [0.3389    , 0.3389    , 0.32889998, ..., 0.67359996,
          0.6591    , 0.6591    ],
         [0.36089998, 0.33879998, 0.33879998, ..., 0.67359996,
          0.56839997, 0.56839997],
         ...,
         [0.3225    , 0.3231    , 0.3172    , ..., 0.66209996,
          0.66209996, 0.5502    ],
         [0.31829998, 0.31829998, 0.3277    , ..., 0.71389997,
          0.71389997, 0.6436    ],
         [0.3306    , 0.3117    , 0.31039998, ..., 0.6608    ,
          0.62409997, 0.6357    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-07-27
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3136    , 0.3051    , 0.3035    , ..., 0.6228    ,
          0.5385    , 0.5051    ],
         [0.3046    , 0.29999998, 0.2953    , ..., 0.617     ,
          0.58349997, 0.58349997],
         [0.3191    , 0.3107    , 0.3107    , ..., 0.63989997,
          0.52459997, 0.52459997],
         ...,
         [0.2589    , 0.2652    , 0.2674    , ..., 0.6885    ,
          0.6885    , 0.5043    ],
         [0.2589    , 0.26479998, 0.26729998, ..., 0.636     ,
          0.53459996, 0.5481    ],
         [0.25509998, 0.25149998, 0.2625    , ..., 0.6588    ,
          0.5395    , 0.55439997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-08-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.2956    , 0.28329998, 0.2852    , ..., 0.6142    ,
          0.5919    , 0.42859998],
         [0.29839998, 0.294     , 0.28309998, ..., 0.624     ,
          0.45549998, 0.45549998],
         [0.3029    , 0.294     , 0.294     , ..., 0.624     ,
          0.5263    , 0.5263    ],
         ...,
         [0.2022    , 0.24059999, 0.24059999, ..., 0.524     ,
          0.524     , 0.5135    ],
         [0.2184    , 0.22809999, 0.2245    , ..., 0.6376    ,
          0.524     , 0.4908    ],
         [0.2279    , 0.21949999, 0.21949999, ..., 0.49789998,
          0.46359998, 0.4754    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2020-08-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.5039    , 0.5137    , 0.49289998, ..., 0.7247    ,
          0.6464    , 0.38729998],
         [0.5564    , 0.5227    , 0.5227    , ..., 0.701     ,
          0.5259    , 0.5259    ],
         [0.5721    , 0.58779997, 0.58779997, ..., 0.4025    ,
          0.33499998, 0.33499998],
         ...,
         [0.5689    , 0.5689    , 0.61509997, ..., 0.79289997,
          0.79289997, 0.518     ],
         [0.5907    , 0.587     , 0.61509997, ..., 0.687     ,
          0.59389997, 0.532     ],
         [0.4201    , 0.4199    , 0.4585    , ..., 0.6235    ,
          0.5595    , 0.49929997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-05-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.4466    , 0.4412    , 0.40559998, ..., 0.466     ,
          0.4675    , 0.38259998],
         [0.45859998, 0.4538    , 0.3983    , ..., 0.466     ,
          0.38619998, 0.38619998],
         [0.47059998, 0.46559998, 0.46559998, ..., 0.3751    ,
          0.41369998, 0.41369998],
         ...,
         [0.3788    , 0.403     , 0.42209998, ..., 0.5765    ,
          0.5765    , 0.48139998],
         [0.4276    , 0.403     , 0.4468    , ..., 0.6254    ,
          0.5765    , 0.547     ],
         [0.41099998, 0.3463    , 0.3811    , ..., 0.6362    ,
          0.5912    , 0.55259997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-06-10
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3761    , 0.35009998, 0.33519998, ..., 0.3572    ,
          0.3319    , 0.3788    ],
         [0.36789998, 0.3506    , 0.3268    , ..., 0.3723    ,
          0.3575    , 0.3575    ],
         [0.37559998, 0.36319998, 0.36319998, ..., 0.3687    ,
          0.3247    , 0.3247    ],
         ...,
         [0.2773    , 0.2782    , 0.28599998, ..., 0.66569996,
          0.66569996, 0.487     ],
         [0.2812    , 0.2706    , 0.28599998, ..., 0.5714    ,
          0.5604    , 0.54359996],
         [0.2818    , 0.2706    , 0.2764    , ..., 0.6457    ,
          0.5749    , 0.5711    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-06-26
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3331    , 0.3205    , 0.3044    , ..., 0.2791    ,
          0.2777    , 0.3906    ],
         [0.3247    , 0.3205    , 0.29909998, ..., 0.3452    ,
          0.31309998, 0.31309998],
         [0.3354    , 0.3247    , 0.3247    , ..., 0.3567    ,
          0.34199998, 0.34199998],
         ...,
         [0.2884    , 0.2884    , 0.2965    , ..., 0.60789996,
          0.60789996, 0.5512    ],
         [0.29839998, 0.2884    , 0.2965    , ..., 0.6081    ,
          0.5819    , 0.5789    ],
         [0.29839998, 0.292     , 0.2986    , ..., 0.5852    ,
          0.5941    , 0.5789    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-07-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.28599998, 0.289     , 0.2713    , ..., 0.274     ,
          0.24      , 0.2791    ],
         [0.282     , 0.289     , 0.2713    , ..., 0.28509998,
          0.2832    , 0.2832    ],
         [0.282     , 0.2785    , 0.2785    , ..., 0.31219998,
          0.27809998, 0.27809998],
         ...,
         [0.24479999, 0.24579999, 0.24579999, ..., 0.5888    ,
          0.5888    , 0.525     ],
         [0.24689999, 0.2403    , 0.2367    , ..., 0.52739996,
          0.49929997, 0.525     ],
         [0.23779999, 0.23779999, 0.2367    , ..., 0.4736    ,
          0.41369998, 0.39999998]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-07-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.26139998, 0.2665    , 0.2633    , ..., 0.2992    ,
          0.2695    , 0.3396    ],
         [0.26139998, 0.2595    , 0.2453    , ..., 0.2914    ,
          0.304     , 0.304     ],
         [0.257     , 0.26389998, 0.26389998, ..., 0.2935    ,
          0.304     , 0.304     ],
         ...,
         [0.2272    , 0.23709999, 0.23779999, ..., 0.6392    ,
          0.6392    , 0.49589998],
         [0.22639999, 0.2358    , 0.2245    , ..., 0.5755    ,
          0.60789996, 0.49589998],
         [0.2278    , 0.2183    , 0.2311    , ..., 0.5076    ,
          0.50699997, 0.4755    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-08-13
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.2483    , 0.2428    , 0.2286    , ..., 0.4034    ,
          0.2827    , 0.3074    ],
         [0.2418    , 0.2366    , 0.23519999, ..., 0.3392    ,
          0.2769    , 0.2769    ],
         [0.2463    , 0.2366    , 0.2366    , ..., 0.2884    ,
          0.2641    , 0.2641    ],
         ...,
         [0.22219999, 0.20369999, 0.21069999, ..., 0.6896    ,
          0.6896    , 0.53889996],
         [0.22219999, 0.20699999, 0.2021    , ..., 0.6547    ,
          0.5714    , 0.51379997],
         [0.21149999, 0.2022    , 0.19999999, ..., 0.6113    ,
          0.5559    , 0.5029    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2021-08-29
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.5281    , 0.5517    , 0.5517    , ..., 0.47849998,
          0.353     , 0.3206    ],
         [0.5956    , 0.5956    , 0.5742    , ..., 0.38189998,
          0.4553    , 0.4553    ],
         [0.6002    , 0.5956    , 0.5956    , ..., 0.38189998,
          0.3353    , 0.3353    ],
         ...,
         [0.4312    , 0.46469998, 0.47509998, ..., 0.5861    ,
          0.5861    , 0.5861    ],
         [0.4312    , 0.4219    , 0.47739998, ..., 0.56409997,
          0.4587    , 0.394     ],
         [0.4606    , 0.44079998, 0.4651    , ..., 0.4684    ,
          0.4684    , 0.4105    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-05-25
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.4754    , 0.47349998, 0.4324    , ..., 0.51739997,
          0.4084    , 0.3692    ],
         [0.5453    , 0.5739    , 0.4576    , ..., 0.41189998,
          0.4084    , 0.4084    ],
         [0.6189    , 0.62009996, 0.62009996, ..., 0.3698    ,
          0.4096    , 0.4096    ],
         ...,
         [0.3021    , 0.34989998, 0.3854    , ..., 0.5815    ,
          0.5815    , 0.5667    ],
         [0.3356    , 0.34989998, 0.3854    , ..., 0.682     ,
          0.6777    , 0.56769997],
         [0.37449998, 0.3356    , 0.359     , ..., 0.6737    ,
          0.6777    , 0.6166    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-06-10
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3599    , 0.3716    , 0.36519998, ..., 0.4274    ,
          0.35349998, 0.4064    ],
         [0.39569998, 0.3872    , 0.3793    , ..., 0.401     ,
          0.3752    , 0.3752    ],
         [0.3994    , 0.4121    , 0.4121    , ..., 0.40879998,
          0.4321    , 0.4321    ],
         ...,
         [0.29999998, 0.3109    , 0.32189998, ..., 0.6532    ,
          0.6532    , 0.5092    ],
         [0.31309998, 0.2972    , 0.3191    , ..., 0.54469997,
          0.54469997, 0.5108    ],
         [0.3277    , 0.3083    , 0.31599998, ..., 0.51309997,
          0.5282    , 0.55619997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-06-26
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.30519998, 0.30159998, 0.2868    , ..., 0.262     ,
          0.2606    , 0.3272    ],
         [0.315     , 0.3177    , 0.2944    , ..., 0.3557    ,
          0.3317    , 0.3317    ],
         [0.3329    , 0.3294    , 0.3294    , ..., 0.362     ,
          0.3349    , 0.3349    ],
         ...,
         [0.2316    , 0.2561    , 0.26229998, ..., 0.51519996,
          0.51519996, 0.484     ],
         [0.2543    , 0.24939999, 0.26909998, ..., 0.6109    ,
          0.65059996, 0.5219    ],
         [0.258     , 0.24139999, 0.26119998, ..., 0.5164    ,
          0.4998    , 0.5172    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-07-12
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.2627    , 0.25939998, 0.2553    , ..., 0.2541    ,
          0.2443    , 0.3182    ],
         [0.26119998, 0.2649    , 0.2529    , ..., 0.3214    ,
          0.3216    , 0.3216    ],
         [0.2683    , 0.2713    , 0.2713    , ..., 0.31329998,
          0.2935    , 0.2935    ],
         ...,
         [0.26299998, 0.2683    , 0.2683    , ..., 0.61289996,
          0.61289996, 0.5198    ],
         [0.2809    , 0.26139998, 0.27719998, ..., 0.601     ,
          0.4712    , 0.5198    ],
         [0.276     , 0.26139998, 0.27719998, ..., 0.47309998,
          0.41959998, 0.459     ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-07-28
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.2434    , 0.2579    , 0.2579    , ..., 0.26839998,
          0.2297    , 0.31509998],
         [0.2528    , 0.26209998, 0.25669998, ..., 0.26839998,
          0.3085    , 0.3085    ],
         [0.2681    , 0.2681    , 0.2681    , ..., 0.3351    ,
          0.3347    , 0.3347    ],
         ...,
         [0.2135    , 0.2234    , 0.2208    , ..., 0.5748    ,
          0.5748    , 0.5306    ],
         [0.2198    , 0.21519999, 0.2166    , ..., 0.578     ,
          0.5748    , 0.5276    ],
         [0.2173    , 0.21599999, 0.2235    , ..., 0.5766    ,
          0.47059998, 0.48079997]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-08-13
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area,
 <xarray.DataArray 'NDVI' (date: 1, y: 390, x: 852)> Size: 1MB
 array([[[0.3922    , 0.4044    , 0.4436    , ..., 0.24569999,
          0.242     , 0.3066    ],
         [0.37809998, 0.3847    , 0.428     , ..., 0.2627    ,
          0.3211    , 0.3211    ],
         [0.37129998, 0.40719998, 0.40719998, ..., 0.3037    ,
          0.2997    , 0.2997    ],
         ...,
         [0.2043    , 0.21159999, 0.21339999, ..., 0.5877    ,
          0.5877    , 0.39769998],
         [0.20639999, 0.20199999, 0.2074    , ..., 0.5326    ,
          0.3971    , 0.3372    ],
         [0.20549999, 0.20709999, 0.2113    , ..., 0.4637    ,
          0.3579    , 0.3372    ]]], shape=(1, 390, 852), dtype=float32)
 Coordinates:
     band         int64 8B 1
   * x            (x) float64 7kB -103.0 -103.0 -103.0 ... -101.2 -101.2 -101.2
   * y            (y) float64 3kB 43.8 43.79 43.79 43.79 ... 42.99 42.99 42.99
     spatial_ref  int64 8B 0
   * date         (date) datetime64[ns] 8B 2022-08-29
 Attributes:
     units:          NDVI
     AREA_OR_POINT:  Area]
InĀ [9]:
# Combine NDVI images from all dates
da = xr.combine_by_coords(ndvi_das, coords=['date'])
/tmp/ipykernel_2071/2757418632.py:2: FutureWarning: In a future version of xarray the default value for compat will change from compat='no_conflicts' to compat='override'. This is likely to lead to different results when combining overlapping variables with the same name. To opt in to new defaults and get rid of these warnings now use `set_options(use_new_combine_kwarg_defaults=True) or set compat explicitly.
  da = xr.combine_by_coords(ndvi_das, coords=['date'])
/tmp/ipykernel_2071/2757418632.py:2: FutureWarning: In a future version of xarray the default value for compat will change from compat='no_conflicts' to compat='override'. This is likely to lead to different results when combining overlapping variables with the same name. To opt in to new defaults and get rid of these warnings now use `set_options(use_new_combine_kwarg_defaults=True) or set compat explicitly.
  da = xr.combine_by_coords(ndvi_das, coords=['date'])
InĀ [35]:
#calculate difference in NDVI
ndvi_diff = (
    da.sel(date=slice('2018')).mean('date').NDVI

    - da.sel(date=slice('2016')).mean('date').NDVI
)

#Plot the difference in NDVI
combined_plot1 =(ndvi_diff.hvplot(x='x', y='y', cmap='PiYG',
                  geo=True, 
                  title='Change in NDVI @ Pine Ridge Indian Community\n2018 vs. 2016')
*
prir_gdf.hvplot(geo=True, fill_color=None, line_color='black'))

combined_plot1
Out[35]:
InĀ [36]:
#calculate difference in NDVI
ndvi_diff2 = (
    da.sel(date=slice('2022')).mean('date').NDVI

    - da.sel(date=slice('2019')).mean('date').NDVI
)

#Plot the difference in NDVI
combined_plot2 =(ndvi_diff2.hvplot(x='x', y='y', cmap='PiYG',
                  geo=True, 
                  title='Change in NDVI @ Pine Ridge Indian Community\n2022 vs. 2019')
*
prir_gdf.hvplot(geo=True, fill_color=None, line_color='black'))

combined_plot2
Out[36]:
InĀ [Ā ]:
# Compute mean annual July NDVI for 2016-2018
ndvi_july_1618 = (
                    da.sel(date=slice('2016','2018'))
                    .where(da.date.dt.month == 7, drop=True)
                    .groupby('date.year')
                    .mean('date')
                    .NDVI)

print(ndvi_july_1618.head())

# Compute mean annual July NDVI for 2019-2022
ndvi_july_1922 = (
                    da.sel(date=slice('2019','2022'))
                    .where(da.date.dt.month == 7, drop=True)
                    .groupby('date.year')
                    .mean('date')
                    .NDVI)

print(ndvi_july_1922.head())
<xarray.DataArray 'NDVI' (year: 4, y: 5, x: 5)> Size: 400B
array([[[0.3188    , 0.3272    , 0.3464    , 0.35399997, 0.35399997],
        [0.3189    , 0.31755   , 0.31904998, 0.31904998, 0.33069998],
        [0.31379998, 0.31835   , 0.31835   , 0.33195   , 0.3423    ],
        [0.3226    , 0.3226    , 0.3252    , 0.3363    , 0.33269998],
        [0.333     , 0.33185   , 0.3401    , 0.34595   , 0.34595   ]],

       [[0.30275   , 0.30285   , 0.3083    , 0.30885   , 0.30885   ],
        [0.30065   , 0.30525   , 0.3083    , 0.3083    , 0.30885   ],
        [0.3172    , 0.3164    , 0.3164    , 0.31839997, 0.31695   ],
        [0.3247    , 0.3247    , 0.31629997, 0.31375   , 0.32      ],
        [0.3123    , 0.3213    , 0.31440002, 0.3171    , 0.3171    ]],

       [[0.7323    , 0.7112    , 0.69145   , 0.60405004, 0.60405004],
        [0.7336    , 0.72685003, 0.71844995, 0.71844995, 0.69275   ],
        [0.7271    , 0.71615   , 0.71615   , 0.7053    , 0.7014    ],
        [0.66565   , 0.66565   , 0.68595   , 0.67345   , 0.70005   ],
        [0.64695   , 0.66279995, 0.65489995, 0.65744996, 0.65744996]],

       [[0.70585   , 0.7014    , 0.70645   , 0.67149997, 0.67149997],
        [0.7084    , 0.70079994, 0.70419997, 0.70419997, 0.69365   ],
        [0.69200003, 0.70079994, 0.70079994, 0.68964994, 0.70489997],
        [0.70554996, 0.70554996, 0.70554996, 0.70869994, 0.69914997],
        [0.66335   , 0.6727    , 0.67695   , 0.68104994, 0.68104994]]],
      dtype=float32)
Coordinates:
  * x            (x) float64 40B -103.0 -103.0 -103.0 -103.0 -103.0
  * y            (y) float64 40B 43.8 43.79 43.79 43.79 43.79
    band         int64 8B 1
    spatial_ref  int64 8B 0
  * year         (year) int64 32B 2016 2017 2018 2019
Attributes:
    units:          NDVI
    AREA_OR_POINT:  Area
<xarray.DataArray 'NDVI' (year: 4, y: 5, x: 5)> Size: 400B
array([[[0.70585   , 0.7014    , 0.70645   , 0.67149997, 0.67149997],
        [0.7084    , 0.70079994, 0.70419997, 0.70419997, 0.69365   ],
        [0.69200003, 0.70079994, 0.70079994, 0.68964994, 0.70489997],
        [0.70554996, 0.70554996, 0.70554996, 0.70869994, 0.69914997],
        [0.66335   , 0.6727    , 0.67695   , 0.68104994, 0.68104994]],

       [[0.35635   , 0.351     , 0.34995   , 0.3505    , 0.3505    ],
        [0.3467    , 0.3467    , 0.33765   , 0.33765   , 0.33375   ],
        [0.36325   , 0.35219997, 0.35219997, 0.3366    , 0.35525   ],
        [0.3645    , 0.3645    , 0.3466    , 0.3486    , 0.3886    ],
        [0.3777    , 0.37734997, 0.37589997, 0.37510002, 0.37510002]],

       [[0.30955   , 0.30475   , 0.28785   , 0.27225   , 0.27225   ],
        [0.30335   , 0.30475   , 0.2852    , 0.2852    , 0.27749997],
        [0.3087    , 0.30159998, 0.30159998, 0.29764998, 0.29145   ],
        [0.3247    , 0.3247    , 0.3226    , 0.31075   , 0.29685   ],
        [0.32755   , 0.3245    , 0.32549998, 0.33139998, 0.33139998]],

       [[0.28394997, 0.2805    , 0.27104998, 0.26424998, 0.26424998],
        [0.2881    , 0.2913    , 0.27365   , 0.27365   , 0.2644    ],
        [0.3006    , 0.30035   , 0.30035   , 0.27365   , 0.28245   ],
        [0.3094    , 0.3094    , 0.31794998, 0.29904997, 0.3017    ],
        [0.31855   , 0.3186    , 0.31849998, 0.31054997, 0.31054997]]],
      dtype=float32)
Coordinates:
  * x            (x) float64 40B -103.0 -103.0 -103.0 -103.0 -103.0
  * y            (y) float64 40B 43.8 43.79 43.79 43.79 43.79
    band         int64 8B 1
    spatial_ref  int64 8B 0
  * year         (year) int64 32B 2019 2020 2021 2022
Attributes:
    units:          NDVI
    AREA_OR_POINT:  Area
InĀ [48]:
mean_1618 = ndvi_july_1618.mean('year')
mean_1922 = ndvi_july_1922.mean('year')

july_diff = mean_1922 - mean_1618

july_diff.hvplot(
    x='x',
    y='y',
    cmap='PiYG',
    clim=(-0.3, 0.3),
    geo=True,
    title='Change in July NDVI: 2019–2022 minus 2016–2018'
)
Out[48]:
InĀ [44]:
# Plot means between the two time periods

# Compute spatial mean and drop extra columns
ndvi_july_1618_mean = (ndvi_july_1618
                       .mean(['x','y']).to_dataframe()
                       .drop(columns=['band','spatial_ref'], errors='ignore'))


ndvi_july_1922_mean = (ndvi_july_1922
                       .mean(['x','y']).to_dataframe()
                       .drop(columns=['band','spatial_ref'], errors='ignore'))

# Join datasets 
july_ndvi_combined = (
      ndvi_july_1618_mean
            .rename(columns={'NDVI':'NDVI 2016-2018'})
            .join(ndvi_july_1922_mean
                  .rename(columns={'NDVI':'NDVI 2019-2022'}),
                  how='outer')
            .reset_index())

# Plot as lines
july_ndvi_combined.hvplot(
      x='year',
      y=['NDVI 2016-2018','NDVI 2019-2022'],
      title='Mean July NDVI 2016-18 vs 2019-22 in PRIC',
      ylabel='Mean July NDVI')
Out[44]: